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Noneqvulibrium  carrier  distributions  in  microcavity  lasers  are  computed  by  solution  of  a  quantum  Boltzmann 
equation  that  includes  carrier-earner,  camer-phonon,  and  camer-photon  scattering  as  well  as  the  pump 
process.  A  significant  heating  of  the  earner  plasma  is  observed  as  a  consequence  of  the  Pauli  blocking  of  earner 
injection  and  the  removal  of  cold  earners  through  the  process  of  stimulated  recombination. 


The  emission  characteristics  of  semiconductor  lasers 
are  governed  by  the  frequency-dependent  modal  gain 
and  spontaneous  emission.  In  microscopic  theory, 
gain  and  spontaneous  emission  are  explicit  functions 
of  the  momentum  probability  distributions  of  the 
electrons  and  holes.1  These  carrier  distributions 
are  determined  by  the  interplay  between  the  pump 
process,  intraband  carrier  scattering,  and  stimulated 
and  spontaneous  recombination. 

If  the  carrier  generation  occurs  through  carrier 
injection,  the  carriers  can  be  assumed  to  have  ther- 
malized  with  the  lattice  by  the  time  they  arrive  at 
the  active  region.  Hence  the  pump-generated  dis¬ 
tribution  function  can  be  approximated  by  a  quasi¬ 
equilibrium  Fermi-Dirac  distribution  at  the  lattice 
temperature.  At  the  beginning  of  the  pump  process, 
all  carrier  states  in  the  previously  unexcited  laser 
region  are  empty  and  can  be  populated.  However, 
as  soon  as  the  bands  are  partially  filled,  Pauli  block¬ 
ing  becomes  effective,  and  the  pumping  is  gradually 
blocked.  Since  this  blocking  is  more  efficient  for 
energetically  lower  states,  the  main  kinetic  energy 
of  excited  carriers  increases. 

During  the  pump  process  the  carriers  in  the  laser 
region  are  constantly  redistributed  into  thermal  dis¬ 
tributions  by  the  very  fast  carrier-carrier  Coulomb 
scattering  that  leads  to  scattering  times  of  -100  fs.2 
However,  since  this  elastic-scattering  process  does 
not  dissipate  kinetic  energy,  the  carrier  system 
equilibrates  at  a  temperature  that  is  higher  than 
that  of  the  lattice.  Energy  from  the  carrier  system  is 
transferred  to  the  lattice  through  the  carrier -phonon 
coupling.  Even  though  electron -hole  scattering 
tends  to  equalize  the  temperature  between  these 
two  subsystems,  the  different  pump  blocking  may 
result  in  different  effective  temperatures  for  electrons 
and  holes. 

Stimulated  recombination  of  carriers  that  is  due 
to  the  lasing  process  occurs  in  the  spectral  region 
in  which  the  carrier  system  is  inverted.  Since  this 
inversion  region  is  located  just  above  the  (renor¬ 
malized)  band  gap,  the  stimulated  emission  removes 
electrons  and  holes  with  below-average  kinetic  en- 
ergy,  which  further  increases  the  carrier  tempera¬ 
ture.  Well  above  threshold,  when  the  stimulated 
recombination  time  is  comparable  to  or  shorter  than 
the  intraband  carrier  scattering,  significant  devia¬ 


tions  of  the  actual  carrier  distributions  from  the 
quasi-equilibrium  Fermi-Dirac  statistics  have  been 
found.3 

In  this  Letter  we  investigate  nonequilibrium  carr¬ 
ier  distributions  in  semiconductor  lasers.  Our  anal¬ 
ysis  is  performed  for  microcavity  lasers,  since  these 
permit  only  single-longitudinal-mode  operation.4 
However,  the  conclusions  are  also  valid  for  conven¬ 
tional  long-cavity  lasers.  We  study  the  temperature 
changes  of  the  carrier  system  and  their  influence 
on  the  stationary  and  dynam  -al  laser  operation  by 
solving  a  Boltzmann  equation  for  the  wave-vector- 
dependent  occupation  functions  of  carriers  fl'k: 

j-t  ft  -  (s*a)«« t  +  (s*°). *m  +  +  (s;)»»»p .  (i) 

The  different  rates  on  the  right-hand  side  of  this 
equation  describe  the  carrier  scattering  (electron- 
electron,  electron-hole,  hole-hole,  carrier -LO- 
phonon,  carrier- LA-phonon),  the  stimulated  and 
spontaneous  recombination  of  carriers  that  is  accom¬ 
panied  by  photon  emission,  as  well  as  the  pump 
process.3  The  pump  process  is  modeled  by 


(Si)  P 


jy  Pj(i  -  fj) 
ed  ' 


(2) 


s 


where  j,  e,  d,  and  -q  are  the  pump  current,  electron 
charge,  thickness  of  the  active  region,  and  geometric 
quantum  efficiency,  respectively.  The  Fermi-Dirac 
function  P*  «  (1  +  exp[(l/*D(«J  -  /iJI}'1  describes 
the  energy  distribution  of  carriers  arriving  in  the 
active  region.  The  chemical  potential  of  electrons 
and  holes  depends  on  the  applied  electric  field, 
eih  are  the  free- particle  energies  of  electrons  and 
holes,  and  T  is  the  lattice  temperature. 

Without  Pauli  blocking,  (1  -  /*)  =*  1,  and  the  total 
pump  rate  1/V  Xj,(Si)pamg  reduces  to  jy/ed.  The 
Pauii-blocking  factor  decreases  the  pump  efficiency 
and  modifies  the  energy  distribution  of  excited  car¬ 
riers.  It  is  therefore  meaningful  to  introduce  an 
effective  quantum  efficiency  for  electrons  and  holes, 
which  is  the  product  of  the  geometric  quantum  effi¬ 
ciency  and  the  Pauii-blocking  factor  (yi)tg  —  y(l  - 
ft).  This  A -dependent  quantum  efficiency  is  dose  to 
zero  for  the  electrons  in  the  spectral  region  of  the 
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Fig.  1.  Time-resolved  evolution  of  carrier  density  (thick 
curves),  laser  intensity  (thin  curves),  and  pump  rate 
normalized  to  j  rj/ed  for  various  energy  distribution 
functions  of  carriers  arriving  in  the  active  region 
Pa  «  {1  +  expt(l/AT)(e*  -  Mo)]}"1  with  lattice  tempera¬ 
ture  T  -  300  K  and  density  NP  -  1/V  £*  Pk.  Zen  time 
corresponds  to  the  onset  of  the  pump  process. 

gain  and  approaches  the  usual  quantum  efficiency 
for  higher-lying  k  values. 

To  compute  the  nonequilibrium  carrier  distribu¬ 
tions  under  lasing  conditions  we  have  to  solve  Eq.  (1) 
together  with  the  equation  for  the  spectral  light  in¬ 
tensity  /(to)  *  l/V  £,/,(«), 

yt  /,(«)  -  £  g(to)  -  j-  /,(<*)  +  W(to)S,(to). 


where  optical  gain  g{o>)  and  spectral  rate  of 
spontaneous  emission  W(o>)  depend  directly  on  the 
carrier-momentum-dependent  occupation  functions 
of  electrons  and  holes.3  The  cavity  structure  of  the 
laser  device  modifies  the  density  of  state  of  the 
various  light  modes.5  In  Eq.  (3)  this  leads  to  a 
modified  loss  rate  */Ff  and  medium  refractive  index 
n(to)/F,,  where  F,  describes  the  eigenmodes  of  the 
laser  cavity.  The  line  shape  function  S,  (<w )  describes 
the  connection  between  photon  wave  vectors  q  and 
light  frequencies  to  according  to  the  modified  photon- 
dispersion  relation,  where  l/V  £,  S,(w)  yields  the 
photon  density  of  state. 

We  numerically  solved  the  coupled  Eqs.  (1)  and 
(3)  for  different  pump  situations,  where  the  initial 
condition  is  the  unexcited  semiconductor  and  zero  in¬ 
tensity  inside  the  cavity.  We  assume  that  the  pump 
rate  is  ramped  up  to  its  stationary  value  within 
the  first  picosecond.  In  Figs.  1  and  2  we  compare 
the  cases  of  different  mean  kinetic  energies  of  in¬ 
jected  carriers,  i.e.,  various  chemical  potentials 
in  P*.  for  otherwise  the  same  conditions  to  distin¬ 
guish  clearly  the  influence  of  pump  blocking  from  the 


fact  that  for  the  same  blocking  an  increase  of  the 
pump  current  would  also  increase  the  temperature. 
Figure  1  shows  the  time-resolved  density,  laser  in¬ 
tensity,  and  pump  rate  for  three  different  sets  of 
M.  *  according  to  the  given  densities  .VP  -  1/  V  Pk . 
The  corresponding  occupation  functions  of  electrons 
are  displayed  in  Fig.  2  for  the  stationary  situation. 
For  higher  chemical  potential  [higher  Np-  Fig.  2(ai], 
the  effective  pump  function  ( 1  -  )Pk  has  a  peak  at 
larger  wave  numbers,  indicating  the  preferred  pump¬ 
ing  of  energetically  higher  states.  Correspondingly, 
we  find  a  strong  deviation  of  the  actual  carrier  oc¬ 
cupation  fk  from  a  lattice-temperature  Fermi -Dirac 
distribution  with  the  same  carrier  density.  For  low 
chemical  potentials  [NP  *3  x  1017  cm'3;  Fig.  2(c)], 
the  pumping  of  carriers  according  to  P*  is  restricted 
to  the  region  of  effective  blocking.  This  results  in  a 
lower  carrier  density  (Fig.  1,  top)  because  of  the  less 
effective  pumping  (Fig.  1,  bottom).  Furthermore,  it 
yields  a  lower  effective  carrier  temperature  owing  to 
preferred  pumping  of  energetically  lower  states.  In 
the  presence  of  the  running  laser  mode,  the  effective 
carrier  temperature  increases  further  as  a  conse¬ 
quence  of  stimulated  recombination  of  cold  carriers  in 
a  wave-vector  region  around  the  laser  mode  [marked 
by  arrows  in  Figs.  2(a)-2(c)]. 

In  Figs.  3  and  4  we  now  analyze  the  situation 
in  which  both  the  pump  current  and  the  chemi- 
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Fig.  2.  Momentum-dependent  electron  distribution  Auc¬ 
tions  (solid  curves)  for  the  conditions  of  Fig.  1  after 
steady-state  conditions  have  been  reached.  The  pump 
densities  are  (a)  NP  -  10“  cm*5,  (b)  Np  *  3  x  10“  cm  , 
(c)  Np  »  3  x  10“  cm'3.  The  dotted  curve  shows  Pa,  the 
dotted-dashed  curve  is  Pt(l  -  /*),  and  the  short-dotted 
curve  shows  a  Fermi -Dirac  distribution  function  ft  with 
the  same  density  as  fa  and  lattice  temperature  T  -  300  K. 
The  effective  carrier  temperaturse  are  (a)  435  K,  (b) 
390  K,  and  (c)  355  K.  .  . 
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Fig.  3.  Time-resolved  evolution  of  carrier  density  (thick 
curves),  laser  intensity  (thin  curves),  and  pump  rate  nor¬ 
malized  to  j  Tiled  for  various  pump  currents  j  and  carrier 
distributions  P*.  Zero  time  corresponds  to  the  onset  of 
the  pump  process. 
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Fig.  4.  Stationary  electron  distribution  functions  (solid 
curves)  corresponding  to  Fig.  3  for  (s)AJ>  -  3x  10“  cm-3, 
(b)  Np  _  10“  cm'3,  and  (c)  Np  -  3  x  10“  cm'3.  The 
different  curves  are  ss  in  Fig.  2.  The  effective  carrier 
temperatures  are  (a)  570  K,  (b)  440  K,  and  (c)  340  K. 

cal  potential  of  the  energy  distribution  of  excited 
earners  P*  change  simultaneously.  For  this  pur¬ 


pose  we  assume  that  Np  corresponds  to  the  density 
that  would  be  generated  by  the  pump  process  alone, 

i. e.,  the  density  generated  in  the  absence  of  optical 
transitions.  Assuming  for  simplicity  a  nonradiative 
lifetime  ?or  of  a  process  proportional  to  the  carrier 
density,  we  see  that  Np  follows  from  the  balance 
j  v/ed  —  Np/rnr.  Note  that  the  additional  consider¬ 
ation  of  other  nonradiative  processes,  such  as  Auger 
recombination,  leads  to  only  a  shift  of  the  operational 
point  of  the  laser.  Figure  3  shows  carrier  density, 
laser  intensity,  and  pump  rate  for  various  pump 
currents  j  riled  *  8.7  x  1015,  2.9  x  1018,  and  8.7  x 
10“  cm'3  ps'1  and  the  corresponding  pump  densities 
NP  -  3  x  1018,  1  x  10“,  and  3  x  10“  cm'3.  Since 
we  now  identify  higher  Np  and,  correspondingly, 
lower  pump  blocking  with  increased  pump  current 

j,  the  changes  of  the  stationary  laser  intensity  are 
even  more  pronounced  than  in  Fig.  1.  In  Fig.  4  we 
see  that  the  preferred  pumping  of  the  energetically 
higher  states,  accompanied  by  a  higher  total  pump 
rate  that  increases  the  effective  carrier  that  tem¬ 
perature  even  more  than  in  the  case  of  Fig.  2.  It 
should  be  noted  that  the  gain  changes  that  occur  as 
a  consequence  of  carrier  heating  and  hole  burning 
prevent  a  clamping  of  the  total  carrier  density  for 
above-threshold  pumping  levels.  This  indicates  that 
the  carrier  density  is  not  a  good  measure  for  the 
modal  gain  under  the  discussed  conditions.3 
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